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The  Radon  transform  is  well-known  as  the  basis  To f"  medical 
computed  tomography,  which  produces  a  two-dimensional  (2-D)  map  of  x-ray 
attenuation  coefficient  within  the  body  from  i-D  projection  data.  The 
general  Radon  transform  operates  on  a  data  set  of  M-dimensions ,  reducing 
the  data  to  a  sequence  of  1-U  data  via  integration  over  (M-1)  dimensions. 
The  existence  of  very  efficient  1-D  signal  processing  components,  such  as 
surface  acoustic  wave  (SAW)  filters  or  charge-coupled  devices  (CCDs), 
makes  the  reduction  in  dimensionality  very  useful  for  many  processing 
applications.  Many  operations  can  be  performed  rapidly  on  2-D  or  3-D 
data  by  performing  the  Radon  transform  and  processing  in  1-0  with  these 
sophisticated  processors.  Such  operations  as  Fourier  transformation, 
convolution,  correlation,  bandwidth  compression,  and  computation  of  image 
moments,  ambiguity  function:  or  Wlgner  distribution  functions  can  be 
carried  out  in  this  manner. 


This  document  is  a  report  of  technical  progress  made  during  tne 
second  year  of  a  proposed  three-year  program  to  Investigate  the 
application  of  the  Radon  transform  to  signal  processing.  The  stated 
goals  of  the  program  are: 

(1)  theoretical  investigation  of  the  role  of  the  Radon  transform  in 
signal  processing; 

(2)  construction  of  a  practical  system  for  2-D  spectral  analysis  and 
image  filtering  at  races  of  at  least  3U  frames  per  second; 

(J)  proof-of-principie  experiments  for  other  processing  operations; 

(‘♦)  determination  of  the  feasibility  of  3-D  Radon  processing. 
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In  Che  second  year,  the  specific  research  casks  co  be  performed 

were: 

(1)  concinue  cheorecical  analysis  of  che  Radon  cransform  and  ics 
applicaCions  in  signal  processing; 

(2)  finish  conscruccion  and  cescing  of  an  opcical  Radon  cransformer 
based  on  a  flying-line  scanner; 

(3)  complece  conscruccion  and  cescing  of  a  SAW  Fourier  cransformer 
including  demonscracion  of  che  abilicy  Co  compuce  complex 
cransforms; 

(A)  use  che  opcical  Eladon  cransformer  and  SAW  Fourier  cransformer  for 
2-D  speccrum  analysis; 

(5)  design  and  conscrucc  an  opcical  back-projeccor; 

^6)  design  a  SAW  filcer  wich  a  |v  |-filcer; 

(7)  begin  conscruccion  and  cescing  of  che  |v  |-filcer; 

(3)  invescigace  che  feasibilicy  of  3-D  daCa-processing  syscems  based  on 
speccral  hole  burning; 

(9)  scudy  che  performance  of  a  Radon  correlacion  cracker  by  compucer 
simulacion; 

(lU)  demonscrace  che  principle  of  che  correlacion  cracker  wich  a 


breadboard  experimenc. 


2.  STATUS  OF  THE  RESEARCH 


Progress  made  coward  accomplishing  che  specific  goals  for  the 
second  year  is  as  follows: 

2.1  Theorecical  invesclgacion  of  che  Radon  cransform  applied  Co  signal 
processing. 

The  work  in  chis  area  has  concencraced  on  che  Wigner  discribucion 
funccion  (WDF).  The  WDF,  a  simulcaneous  coordinace-frequency 
represencacion  of  a  funccion,  has  applicacions  in  syscem  descripcion  and 
paccern  recognicion.  Because  ic  is  compucacionally  demanding,  ic  is  noc 
ofcen  implemenced  digically,  especially  for  2-D  inpuc  daca.  Several 
syscems  have  been  proposed  co  compuce  che  WDF  opcically:  eicher  che 
full  2-0  funccion  of  1-D  inpuC  daCa  or  2-D  slices  of  che  4-D  funccion  of 
2-D  inpuc  daca.  We  have  shown  chac  che  Radon  cransform  is  an  effeccive 
cool  for  calculacing  che  WDF  because  of  che  reduccion  co  a  sequence  of 
1-D  Fourier  cransforms  of  an  easily  generaced  funccion.  The  oucpuc  daca 
can  cake  che  form  of  eicher  1-D  lines  or  2-D  slices  of  che  4-D  WDF.  We 
have  successfully  implemenced  chis  operacion  and  published  che  resulcs 
In  Opcical  Engineering.  Invesclgacion  and  implemencacion  of  ocher 
operacions  in  Radon  space,  such  as  2-D  scale  invarianc  paccern 
recognicion  and  2-D  cepscrum  analysis,  is  concini:ing. 


2.2  Compiecion  of  the  optical  Radon  transformer. 

All  of  the  experimental  worR  completed  thus  far  has  been 
accomplished  using  one  version  of  the  fl)ring— line  scanner,  shown 
schematically  in  Fig.  1.  It  produced  a  line  of  he-Ne  laser  light  with  an 
anamorphic  optical  system  and  projected  the  line  onto  the  input  data 
plane.  An  acousto-optic  dragg  cell  driven  by  a  voltage-controlled 
oscillator  enabled  the  line  to  be  scanned  across  the  data  plane, 
resolving  approximately  150  lines  in  an  aperture  of  30  mm  diameter.  A 
single  scan  was  performed  in  10^  s,  yielding  one  line  through  the  Radon 
transform  of  the  input  data.  The  azimuth  of  the  scan  angle  could  be 
varied  with  an  image-rotation  prism,  allowing  the  complete  Radon 
transform  to  be  generated.  The  time  required  to  produce  the  transform 
is  limited  by  the  rotation  speed  of  the  prism  to  approximately  0.3  Hz, 
but  a  new  prism  system  now  under  construction  will  achieve  the  required 
7.5-Hz  rotation  rate,  giving  a  30-Hz  frame  rate.  The  first  version  of 
the  optical  Radon  transformer  wasted  much  light,  requiring  a 
photomultiplier  detector.  A  new,  efficient,  and  compact  system  has  been 
built,  allowing  the  use  of  a  fast  P-I-N  photodiode  detector  and  making 
computation  of  complex  Fourier  transforms  more  convenient. 

2.3.  Completion  of  the  SAW  chirp  Fourier  transformer. 

Production  of  1-D  power  spectra  with  a  SAW  Fourier  transformer  was 
easily  done,  but  implementation  of  the  complete  l-D  complex  Fourier 
transform  was  more  difficult  because  of  the  accurate  timing  required. 
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Flying-Line  Scanner.  The  laser  source  is  imaged  to  a  Line  on  the 
input  plane  by  the  cylindrical  lens.  The  Bragg  cell  is  driven  by  a 
chirp  signal  from  a  voltage-controlled  oscillator,  producing  a  phase 
grating  of  linearly  varying  frequency  that  scans  the  line  across  the 
input.  The  light  transmitted  or  reflected  from  the  input  is 
collected  by  the  detector.  One  scan  across  the  input  produces  one 
projection  of  the  Radon  transform.  Rotation  of  the  scan  azimuth  by 
the  prism  allows  the  production  of  the  complete  Radon  transform. 
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This  problem  has  been  solved,  and  we  are  now  able  co  generate  the 
complex  Fourier  transform  of  real  input  data.  Analog  signals  that  are 
lOu  s  long  can  be  transformed  in  less  than  30 u  s,  yielding  the  real  and 
imaginary  parts  of  the  Fourier  transform  (i.e.,  sine  and  cosine  transform 
respectively)  in  separate  channels.  A  schematic  diagram  of  the  complex 
SAW  chirp  Fourier  transformer  is  shown  in  Fig.  2.  The  resolution  of  the 
transform  is  limited  to  50  data  points  by  the  time-bandwidth  product  of 
one  of  the  SAW  chirp  filters.  Results  of  the  transformation  of  a 
sequence  of  1-D  scans  are  shown  in  Fig.  3,  along  with  a  computer 
simulation.  Performance  of  the  transformer  is  seen  to  be  excellent. 

2.4.  Generation  of  2-D  spectra  with  the  system. 

Using  the  optical  P.adon  transformer  and  SA.W  chirp  Fourier 
transformer,  2-D  power  spectra  have  been  produced  for  several  objects, 
both  transmissive  and  reflective.  Selected  results  are  presented  in 
Fig.  4.  In  addition,  the  technique  has  been  applied  to  a  simulation  of 
Labeyrie  stellar  speckle  interferometry,  by  averaging  a  sequence  of 
simulated  stellar  speckle  patterns  to  produce  diffraction-limited 
information  (i.e.,  the  autocorrelation)  of  the  original  source.  This  wortc 
has  also  been  published  in  Optical  Engineering. 

2.5.  Construction  of  the  optical  bacx-projector. 

The  optical  back-projector  is  composed  of  two  distinct  parts:  the 
high-pass  filter  and  the  optical  projection  system.  The  optical 
projector  is  simple,  straightforward,  ana  has  been  previously 
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Complex  SAW  Chirp  Fourier  Transformer.  The  major  components  of  the 
SAW  transformer  are  the  three  SAW  chirp  filters.  The  impulse 
response  of  such  a  filter  is  h(t)  ■  Rert(t/T)  cos  (ii),t  t  at*). 
Taking  the  plus  sign  in  the  cosine  term  yields  an  upchirp  (output 
frequency  increases  with  time),  while  the  chirp  with  the  negative 
sign  is  a  downchirp.  To  produce  the  complex  Fourier  transform,  the 
input  signal  is  multiplied  by  a  downchirp  in  an  rf  mixer.  The 
product  is  applied  to  a  filter  with  an  upchirp  impulse  response. 
Square-law  detection  of  the  output  signal  yields  the  squared- 
modulus  of  the  Fourier  transform.  Post-multiplication  of  the 
filtered  signal  with  a  second  upchirp  in  a  phase  comparator  gives, 
after  low-pass  filtering,  the  real  and  imaginary  parts  of  the 
Fourier  transform  (l.e.,  the  cosine  and  sine  transform  respectively). 
Of  course,  the  finite  windows  of  the  chirps  limit  the  extent  of  the 
Fourier  spectrum  that  is  computed  as  well  as  the  resolution 
obtainable. 
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One-dimensional  complex  Fourier  transforms  with  the  SAW  transformer. 
This  is  a  sequence  of  four  signals  and  complex  transforms  from  the 
SAW  transformer,  along  with  computer  simulations.  The  object  was  a 
grating  (75%  transparent,  25%  opaque)  in  a  circular  aperture  of 
20-mm  diameter.  In  each  case,  the  first  trace  is  the  output  of  the 


flying-line  scanner,  i.e.,  one  projection  of  the  Radon  transform.  The 
second  and  third  traces  are  the  cosine  and  sine  transform  outputs, 
respectively,  of  the  SAW  complex  chirp  transformer,  a)  The  grating 
was  centered  in  the  circular  aperture,  giving  a  symmetric  object 
yielding  a  purely  real  transform,  i.e.,  the  sine  transform  vanishes, 
b)  ,  c),  and  d),  the  grating  was  translated  incrementally  in  the 
aperture,  producing  a  nonzero  sine  transform.  Comparison  of  each 
case  with  the  computer  simulations  shows  excellent  agreement. 
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4a.  The  object  was  a  pair  of  identical  circular  apertures,  and  the  cosine 
fringes  across  the  Airy  disk  pattern  are  clearly  seen. 


4b.  The  object  consisted  of  two  identical  gratings  crossed  at  90®.  The 
fundamental  orders  of  each  grating  and  the  cross-terms  are  clearly 
seen. 


4.  Two-dimensional  spectra  obtained  with  Che  Radon-Fourier  transformer. 
Five  examples  of  two-dimensional  spectra  are  shown.  Figures  a)  to 
c)  show  the  power  spectrum  of  the  object  (i.e.,  squared-magnitude  of 
the  Fourier  transform),  and  d)  and  e)  are  complex  transrorms. 


■^c.  The  object  whose  transform  is  shown  in  c)  was  a  computer-generate 
hologram.  The  Impulse  response  of  the  hologram  is  an  annulus,  as 
clearly  shown. 


4.  Two-dimensional  spectra  obtained  with  the  Radon-Fourier  transforms 
Five  examples  of  two-dimensional  spectra  are  shown.  Figures  a)  t< 
c)  show  the  power  spectrum  of  the  object  (i.e.,  squared-magnitude  c 
the  Fourier  transform),  and  d)  and  e)  are  complex  transforms. 


4d.  and  a.  The  positive  parts  of  the  cosine  and  sine  transforms  of  a 

single  circular  aperture  are  shown,  d)  The  object  is  centered 
in  the  flying-line  scanner,  3aelding  a  purely  real  Fourier 
transform.  In  d.l),  the  positive  part  of  the  cosine  transform 
shows  the  Airy  pattern,  while  in  d.2),  the  sine  transform  is 
nearly  zero.  e)  The  aperture  has  been  translated  in  the 
scanner,  yielding  a  linear  phase  term  and  a  nonzero  sine 
transform.  In  e.l),  the  cosine  transform  is  seen  to  be 
symmetric  and  the  fringes  due  to  the  Linear  phase  term  are 
shown.  The  sine  transform  in  e.2)  is  antisymmetric. 


4.  Two-dimensional  spectra  obtained  with  the  Radon-Fourier  transformer. 
Five  examples  of  two-dimensional  spectra  are  shown.  Figures  a)  to 
c)  show  the  power  spectrum  of  the  object  (i.e.,  squared-magnitude  of 
the  Fourier  transform),  and  d)  and  e)  are  complex  transforms. 


demons  traced.  The  high-pass,  or  |v  |-£ilcer,  can  be  implemented  in 
several  ways;  the  most  desirable  is  a  custom-designed  SAW  filter.  This 
approach  will  be  discussed  below.  Another  feasible  method  utilizes  a 
SAW  convolver  with  the  filter  function  stored  in  a  fast  digital  memory. 
This  technique  has  been  tested  with  a  prototype  digital  memory  at  low 
data  races  and  been  found  to  work  satisfactorily.  However,  at  the 
speeds  required  for  video-rate  operation,  the  digital  memory  is  far  coo 
noisy.  Components  for  an  improved  memory  using  emitter-coupled  logic 
(ECL)  are  on  order  and  should  solve  the  problem. 

2.6  and  2.7.  Design  and  Construction  of  SAW  [v  ]-filcer. 

The  basic  design  for  the  SAW  |v  [-filter  required  for  the  optical 
back-projector  has  been  completed.  It  will  consist  of  two  SAW  filters 
on  the  same  substrate.  Both  filters  will  utilize  chirped  transducers, 
oriented  to  obtain  a  non-dispersive  output.  To  obtain  the  required 
transfer  function,  two  filter  outputs  will  be  needea:  a  triangular 
bandpass  will  be  subtracted  from  a  rectangular  bandpass.  A  schematic 
diagram  of  the  filter  is  snown  in  Fig.  3.  The  preliminary  test  masKs  for 
the  photolithography  have  been  produced,  and  the  rirsc  test  filter  is 


being  constructed. 


Schemaclr  diagram  of  Che  surface  acoustic  wave  |v|-filcer.  To 
allow  for  accurate  nulling  at  the  center  frequence,  the  filter  is 
divided  into  two  parts.  The  incoming  signal  is  split,  with  half 
applied  to  a  filter  with  a  triangular  transfer  function  and  half  to 
a  filter  with  a  near-rectangular  passband.  A  hybrid  phase-shifter 
in  one  arm  compensates  for  differencial  phasing  in  Che  SAW  filters. 
The  outputs  are  subtracted  in  a  180*  hybrid.  The  output  transfer 
function  has  the  required  shape. 


2.8.  Investigate  the  feasibility  of  3-D  data-processing  systems  based  on 
spectral  hole  burning. 

The  goal  for  the  second  year  in  terms  of  optical  storage  of  analog 
data  was  to  demonstrate  wavelength-multiplexed  storage  in  allcali-halide 
crystals.  We  have  investigated  wavelength  tuning  capabilities  of  laser 
diodes  and  test-irradiated  a  crystal  of  lithium  fluoride  (LiF)  to  observe 
spectral  hole-burning.  The  optical  and  electronic  systems  are  not  yet 
in  their  final  forms,  but  early  results  are  promising.  The  optical 
system  is  shown  in  Fig.  6a.  Using  the  lamp  source  alone,  we  have 
observed  Che  absorption  line  to  be  used  for  storage,  designated  as  A  in 
Fig.  6b.  Tuning  the  laser  diode  to  a  wavelength  within  the  absorption 
line  and  recording  the  transmission  as  a  function  of  time,  we  have 
observed  the  burning  of  a  hole  as  shown  in  Fig.  6c.  The  increase  in 
transmission  was  about  L0%  in  40  secs  of  exposure.  No  attempts  have 
been  made  to  determine  how  much  of  the  laser  power  is  passing  through 
the  spectrometer  slit,  but  we  expect  to  be  able  to  increase  the  burning 
rate  significantly  by  modifying  the  optical  system  to  ensure  that  all 
the  power  reaches  the  crystal. 

We  have  also  used  the  laser-diode  drive  current  to  scan  tne 
wavelength  through  the  spectrum  surrounding  the  burned  hole.  Because  of 
the  small  depth  of  the  hole  and  the  fairly  insensitive  detector 
electronics,  we  were  unable  to  observe  conclusively  the  presence  of  the 
burnt  hole.  We  are  currently  oeveloplng  improved  detector  circuitry  to 
reduce  the  amplifier  noise,  and  are  using  derivative  spectroscopy,  rather 


chan  less  sensitive  transmission  spectroscopy,  to  observe  the  holes. 
This  should  allow  us  to  easily  read  hole  depths  of  the  magnitude  we  will 
be  writing. 

2.9.  Study  the  performance  of  a  Radon  correlation-cracker  by  computer 
simul cation. 

2.10.  Demonstrate  the  correlation-cracRer  with  a  breadboard  experiment. 

We  have  found  published  results  of  work  very  similar  to  that  which 

we  planned  to  do  (Yao  and  Adrian,  Appl.  Opt.,  16d7,  1984),  and  so  we 
concentrated  efforts  in  another  area  instead.  We  have  determined  that 
it  is  feasible  to  partially  perform  some  operations  during  che  scanning 
operation.  For  example,  differentiation  of  che  input  data  could  be  done 
by  scanning  with  a  Laplacian  operator  and  back-projecting.  The  Laplacian 
operator  could  be  generated  by  using  three  side-by-side  light  sources  in 
che  scanner,  with  che  center  source  twice  as  bright  and  18l/  out-of¬ 
phase  relative  to  che  ocher  two  sources.  Suitable  sources  include  laser 
diodes,  but  the  available  wavelengths  (>7S0  nm)  are  too  long  to  be  used 
with  our  ttragg  cell  scanner.  Light-emitting  diodes  exist  with  suitable 
wavelengths,  but  Che  maximum  modulation  speed  is  not  sufficient  for  our 
purpose.  A  third  possibility  is  to  use  che  He-Ne  laser  in  che  scanner 
now,  with  fiber  optic  delay  lines  to  properly  phase  che  sources.  This 
approach  to  image  differentiation  should  be  less  sensitive  co  noisy  data 
because  of  che  line-integral  of  the  Radon  transform. 
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6.  FIGURE  CAPTIONS 


Flying-Line  Scanner.  The  laser  source  is  imaged  co  a  line  on  Che 
input  plane  by  the  cylindrical  lens.  The  Bragg  cell  is  driven  by  a 
chirp  signal  from  a  voltage-controlled  oscillator,  producing  a  phase 
grating  of  linearly  varying  frequency  that  scans  the  line  across  the 
input.  The  light  transmitted  or  reflected  from  the  input  is 
collected  by  the  detector.  One  scan  across  the  input  produces  one 
projection  of  the  Radon  transform.  Rotation  of  the  scan  azimuth  by 
the  prism  allows  the  production  of  the  complete  Radon  transform. 


Complex  SAW  Chirp  Fourier  Transformer.  The  major  components  of  the 
SAW  transformer  are  the  three  SAW  chirp  filters.  The  impulse 
response  of  such  a  filter  is  h(t)  ■  Rect(t/T)  cos  (u,t  ±  at*). 
Taking  the  plus  sign  in  the  cosine  term  yields  an  upchirp  (output 
frequency  increases  with  time),  while  the  chirp  with  the  negative 
sign  is  a  downchirp.  To  produce  the  complex  Fourier  transform,  the 
input  signal  is  multiplied  by  a  downchirp  in  an  rf  mixer.  The 
product  is  applied  to  a  filter  with  an  upchirp  Impulse  response. 
Square-law  detection  of  the  output  signal  yields  the  squared- 
modulus  of  the  Fourier  transform.  Post-multiplication  of  the 
filtered  signal  with  a  second  upchirp  in  a  phase  comparator  gives, 
after  low-pass  filtering,  the  real  and  imaginary  parts  of  the 
Fourier  transform  (l.e.,  the  cosine  and  sine  transform  respectively). 
Of  course,  the  finite  windows  of  the  chirps  limit  the  extent  of  the 
Fourier  spectrum  that  Is  computed  as  well  as  the  resolution 
obtainable. 


1-D  Complex  Fourier  Transforms  with  the  SAW  Transformer.  This  is  a 
sequence  of  four  outputs  from  the  flying  line  scanner  together  with 
the  outputs  of  the  chirp  Fourier  transformer.  The  left  hand 
figures  are  the  oscilloscope  traces  from  the  scan  and  transformer, 
while  those  on  the  right  are  computer  simulations  to  demonstrate 
the  accuracy  of  the  computation.  The  object  was  a  grating  (751 
transparent,  25*  opaque)  in  a  circular  aperture,  a)  The  output  of 
the  flying  line  scanner,  i.e.,  one  projection  of  the  Radon  transform. 

b)  The  cosine  and  sine  transforms  of  the  projection  with  the  grating 
centered  In  the  aperture,  resulting  in  a  symmetric  object  with  a 
purely  real  transform,  i.e.,  the  sine  transform  vanishes. 

c)  Examples  of  the  cosine  and  sine  transforms  of  the  projection 
with  the  grating  translated  incrementally  in  the  aperture,  adding  an 
antisymmetric  part  to  the  object,  and  resulting  in  a  nonzero  sine 
transform.  Comparison  with  the  computer  simulations  shows 
excellent  agreement. 


2-D  Spectra  obtained  with  the  Radon-Fourier  Transformer.  Five 
examples  of  2-0  spectra  are  shown,  a),  b),  and  c)  are  the  squared 
moduli  of  the  Fourier  transform  of  the  object,  while  d)  and  e)  are 
complex  transforms,  a)  A  pair  of  Identical  circular  apertures.  The 
cosine  fringes  across  the  Airy  disc  pattern  are  clearly  seen,  b) 
The  two  gratings  crossed  at  90*.  showing  the  fundamental  orders  for 
both  gratings,  plus  the  cross-terms.  c)  A  computer-generated 
hologram  whose  impulse  response  was  an  annulus,  d)  and  e)  Examples 
of  a  single  circular  aperture.  In  d)  the  aperture  was  centered  in 
the  flying-line  scanner,  so  that  the  complex  transform  is  purely 
real.  By  translating  the  object  in  the  flying-line  scanner,  a 
linear  phase  is  added  to  the  transform.  The  sine  transform  is  now 
nonzero,  and  both  the  cosine  and  sine  transforms  show  fringes. 

Schematic  diagram  of  the  surface  acoustic  wave  |v|-filter.  To 
allow  for  accurate  nulling  at  the  center  frequence,  the  filter  is 
divided  into  two  parts.  The  incoming  signal  is  split,  with  half 
applied  to  a  filter  with  a  triangular  transfer  function  and  half  to 
a  filter  with  a  near-rectangular  passband.  A  hybrid  phase-shifter 
in  one  arm  compensates  for  differential  phasing  in  the  SAW  filters. 
The  outputs  are  subtracted  in  a  180*  hybrid.  The  output  transfer 
function  has  the  required  shape. 


Optical  layout  for  the  system  to  store  optical  analog  data  via 
spectral  hole  burning. 


Sample  transmission  spectrum  obtained  for  a  line  to  be  used  for 
storage  (X»833  nm). 

Transmission  of  the  absorption  as  a  function  of  time  during  burn. 
The  rise  in  transmission  is  due  to  depletion  of  the  ground  state 
population. 
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